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ON THE SECOND DERIVATIVES OF THE EXTREMAL-INTEGRAL 
FOR THE INTEGRAL /F(t/; y')dt* 



ARNOLD DRESDEN 

Inteoduction 

It is the purpose of this paper to obtain formulas for the second derivatives 
of the extremal-integral arising in the theory of the integral 

(1) fF{y;y')dt, 

where y , y' are symbols for the w-partite numbers ( j/i , • • • , j/n ) and ( y[ , 
• • • ) y'n) , the ' denoting differentiation with respect to the parameter t , 
expressing these derivatives in terms of " normal solutions " of the Jacobi 
system of differential equations, introduced into the literature by Professor 
Bliss, t The matrix notation used in Bliss's paper will be followed here in 
as far as it is feasible to do so. 

The methods used in this paper are analogous to those of an earlier paper ;t 
for this reason only such details as are essential will be developed here. In 
Section 1 some preliminary facts are put together; in Section 2 the desired 
formulas are derived and discussed. 

1. Preliminary developments 
Using the arc length s as a parameter, we denote by 

(2) y = y{s - si; 2/1, y[) 

the solutions of Euler's differential equations for our problem, which satisfy 
the initial conditions 

2/(0; 2/1, 2/1) = 2/1, 2/' (0; 2/1. 2/1) = 2/1. 

where Si, yii= yn, ■ ■ • , ym) and y'ii= fn, •• • , 2/i») are arbitrary con- 
stants of 1 , n, and n dimensions respectively, subject only to the condition 

* Presented to the Society, December, 1914 and September, 1916. 
t Bliss, these Transactions, vol. 17 (1916), p. 195. 
t Dresden, these Transactions, vol. 17 (1916), p. 425. 
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Sj=i2/ii = 1. The existence of such a solution follows from the general 
existence theorems.* Supposing that we have found an extremal passing 
through Ai{ai) and ^2(a2), not containing a pair of conjugate points, we 
can construct an extremal through the points Pi (2/1) and P2 (2/2) taken suf- 
ficiently near to Ai and A^ respectively. We findf 

(3) y = y{s - si;yi, 9{yi,y2)) = y{s), Si^s^Si, 

where 6 {yi, 2/2 ) and 52(2/1, 2/2 ) are determined so as to satisfy the equations 

(4) 2/(«2 - si;yi, 0) = 2/2, 

i 

The functions y{s)X defined by equations (3) satisfy the following initial 
conditions : 

(5) 2/(«i)=2/i, 2/(^2) =2/2, 2/'(«i) = ^(2/1, 2/2). 

If we write 4> for 2/' ( *2 ) , we can obtain a second form for the equation of 
the extremal Pi P2 , viz. : 

(3a) y = y{s - S2;y2,<i>) ^ y{s) , Si ^ s ^ §2; 

the functions y{s) satisfy the conditions 

(5a) ^(si)=2/i> y{si)=y2, y'isi)=e, y' {s^) = <t) , 

while ^ and S\ satisfy the equations 

(4a) y{si - S2;y2,<l>) = yi, Z)^« = 1. 

i 

Along the extremal Pi P2 taken in the forms (3) or (3a), we compute the 
integral (1). The integral obtained in this way, known as the extremal- 
integral connected with the integral (1), is a function of 2/1 and 2/2; we will 
denote it by T (2/1, 2/2): 

r(2/i,2/2)= rF{y{s),y'is))cb= P F {y{s) , y' {s))ds . 

If we use subscript notation for the derivatives of the function T 
{Tii = dT/dyii, etc., Tii,2j = d'^ T/dyudy2j , etc.). 



* Bolza, Vorlesungen vber Varialionsrechnung, pp. 168-188. 

t Bolza, loc. cit., p. 597. 

t Whenever this can be done without sacrificing clearness, we shall write y (s) in place of 
J/ ( s - si; j/i , e ( j/i , j/2) ) , and F (s) (01 F (y (s); y' (s)); y (s) and F (s) will be used 
with analogous meanings. 
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we find* 

(6) Tiy = - Fn+j{yi; 0) = - Fn+j{Si), T2J = F„+,- (2/2; <j>) = Fn+j (S2) . 

The second derivatives of the function F taken along the extremal Pi P2 
will be denoted by the symbols, used by Bolza,t viz., 

d^F d^F d^F 

dyi dy, - ^^ ' dyi dyl " ^" ' dyl dy'u ~ ^" ' 

For our further work, we introduce the following functions of s : 

Pa = dyi{s) /dy\j , aij = dyiis) /dyij ; 

Pii = dyi{s) /dyij , o-.y = dyi{s) /dy^j . 

From differentiation of both sides of equations (5) and (5a) with respect to 
the variables j/iy and 2/2^ , it follows that these functions satisfy the following 
initial conditions: 

(8) at 5 = Si : Pa = 0, Pij = bif, a^ = 5,7, <r,y = 0; 
(8a) at s = S2 : Pa = 0, p'lj = 5,7; ffij = 6,7, a'ij = 0. 

Furthermore we differentiate both sides of equations (4) and (4a) with 
respect to the variables z/y and y-a , and solve the resulting equations for the 
derivatives of dh and 4>h; we find: 

^(^2)^— = -Dhj{si), A(52)v7- = A/.y(s2), 

(9) _ _ 

A(si)^^=AAy(5i), A(si)f|^.= -Z)Ay(si). 

Here A ( s ) and A ( s ) are determinants of order n + 1 given by the formulas 



A(5) = 



2/' P 

e 



Ais) = 



1 2/ P 
iO <f> 



Dhj is obtained from A by replacing the {h + l)th column by (Xij (i = I , 
• • • , n) , 0; Ahj is_the cof actor of pjh in A. Similar meanings are to be attri- 
buted to Dhj and A^y. 

2. The second derivatives of the extremal-integral^ 

It follows from the general theory that the functions pu and aij form for 
every fixed j a solution of Jacobi's equations;! that y' is also a solution of 

* Bolza, loc. cit., p. 309. 
t Loc. cit., p. 621. 

X The deductions in this section have been materially simplified by suggestions from Pro- 
fessor Bliss, to whom, in particular, the simple determinantal form of the final formulas is due. 
§ See, e. g., Hadamard, Legons sur le calcul des variations, p. 25. 
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these equations has been proved by Bliss,* who has moreover shown that 
to every solution u of Jacobi's equations, there corresponds a unique normal 
solution of the form u — "Ky' , characterized by the condition 

{u -\y')y' = 0, Si^sSs2, 

where X is an uniquely determined n-partite constant. 

The normal solutions which correspond to the solutions py and ffj are given 
by the formulas 

Vj = Pj - \ y' , Kj = <!, - Ky y' , 

where Xy = y' pj and kj = y' (Xj . From equations (8) it follows that these 
normal solutions satisfy the following initial conditions: 

(10) at s = si, rjij = , rj<j = 5ij — di dj , 

and 

(11) at « = si, Uj = 5ii - di dj, tij = — y'i («i) ^} — y'l («i) ^i- 

From equations (10) we conclude that the functions rj, 6 form for i = 1 , • • • , 
n a normal solution whose elements vanish at 5 = Si , together with their first 
derivatives; consequently they vanish identically throughout the interval 
(si S2); i- e., the n normal solutions Vij ij = 1 > • • • , n) are connected by the 
linear relation 

(12) Vnei+Vi2 02+ ■•■ +VinOn =0 (i = l,---,n). 

That these normal solutions do not satisfy any other linear relation becomes 
evident when we notice that the coefficients Cj for such a relation would have 
to satisfy, in consequence of the conditions (10), the equations 

J^ Cj{Sij - di dj) = 0, i. e., d = diJ2 cy 6j- 
3 J 

Hence we would have either £y cy ^y = , from which would follow Cj = , 
or the fact that the c; are proportional to the di. Since there is always one 
non-vanishing 6j , there is no loss in generality, if we suppose, for simplicity of 
notation, 0„ =# 0; in this case the solutions ryi, • • • , rjn-i form a matrix of 
n — 1 linearly independent solutions. 

We differentiate now equations (6) and obtain 

Tii.iy = - QjiiSi) - Z)-RiA(Sl)^ , 

A "i/iy 

Tu, 2y = — X) -^iA (*i) 5^ ' 

and similar formulas for r2,-, ly and r2f, 2y. Substituting equations (9) in 
* Bliss, loc. cit., p. 201. 
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these formulas, we can put them into the form : 
Qij Ri 



Til, 1} = — 







y p 
d 



Tli, 2} — 



e 



y p\ 
ey 



where the matrix y' Tij is obtained from the matrix y' p by replacing the jth 
row hyO,Rii, ••• , Rin . 

To reduce the solutions py and aj , occurring in these formulas, to the corre- 
sponding normal solutions, we subtract from each of the columns containing 
the elements of a solution pj or <Tj , the column containing y' multiplied by the 
proper factor Xy or Kj; this will replace in each of the above determinants, 
p and <r by 7j and f respectively. Multiply now the last n columns in each of 
these modified determinants hy 6i, • • • , On respectively, and place the sum of 
the results in the nth column. In view of equations (12) and of the well- 
known fact that Riy' = , we obtain then the following 

Theoeem. The second derivatives of the extremal integral T are given by the 
formulas 

Qjiisi) Rdsi) 

ry(S2) y'{S2) HiS2) 
n-l 

Tii,2j 12Rihisi)HjHis2) ^ \y'is2)His2)\, 

h~l 



Tu, w = — 



^\y'is2)H(s2)\, 



T2i. ij = JLRihis2)Hjk{si) -e- \y'isi)H{si)\, 



■ 2t, 2/ — 



Qli(S2) ^(52) 

ry(sx) y'{s^) Hisi) 



^\y'{s^)His^)\, 



where H is the matrix of the n — 1 linearly independent normal solutions 
Tji, • • • , r]n-i of Jacobi's equations characterized by conditions (10); Hjh is the 
cof actor of rijh in \ y'H \ ; where H and Hjh are similarly formed from the functions 
Tjjh , obtained by normalizing the functions "pih; and where Ri is the (n — 1 )- 
partite number ( Rn , • • • , iii, „_i ) . * 

It still remains to show that Fi,-, 2/ = T2J, u, and that Th, y and T2i, 2j 
are symmetric in i and j, i. e., we have to prove the following relations: 

n-l 

{A) - Q,,. ( si ) • I y'H I ,=,, + Y,Rih{si)-\ y'Zhj \ .=,, 

A=l 

n-l 

= - Qijisi)- \y'H\,^,^ + YjRjhisi)- \y'Zhi\s=,^, 



*=i 

n-l 



{B) 



\y 



'•H'|.=.i • 12Rih{si)-Hjh{s2) +\y'H\,^„ • ^Rjh{s2)-Hih{si) , 



* The choice of the particular matrix H of n — 1 linearly independent normal solutions 
was based upon the assumption e„ 4= 0; if one of the other e'a is different from zero, an evi- 
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and the relation (A) , obtained from (A) by replacing H and Z hy H and Z 
respectively and by permuting the arguments Si and S2. Here \y'Zhj\ is 
obtained from | y'H \ by replacing the functions j?»a by f jy ( i = 1 , • • • , n); 
and 1 y'Zhj I is similarly formed from | y'H \ . 

To this end we begin by noticing that for any two solutions of Jacobi's 
equations, the relation 

(13) u[Qi) + Rv'] - v[Qu + Ru'] = constant 

holds.* Applying this formula in turn for u = ^i, v = tjj , and for m = f», 

V = ^j, and making use of equations (10) and (11), we obtain the results 

^i I Qvj + Rv'j ]-vjl Qti + R^'i] = Rii ( si ) , 

and 

ry [ Qr* + Rti ] - ra Qh + i?r; i = Qh ( ^i ) - Qh ( ^i ) , 

by means of which formula {A) reduces to 

\y'H\-T. Ui UdQik - Qu] + Z Ui hi [R,k m - Roi m]\ = 0; 

*' *'» |s = «2 

where H'i represents the determinant obtained from | y'll \ by replacing the 

elements of the A;th row by the derivatives of the corresponding elements of 

the gtYi row. This last formula follows however immediately from equation 

(13) if we put u = ^k and v = ^i, where ^,„ (s) , i = 1 , • • • , ra, denotes the 

solution of Jacobi's equations obtained by replacing the gt\i row of | y'H \ ,=sj 

by the ith row of \y'H\ and dividing the result by the former determinant; 

we have then Hgis^) = \y'H\,=,^ • ^'i^jisi). Formula (A) is proved in the 

same way, by the use of the functions ^y, fy, and $y in equation (13), while 

formula (B) follows also from this same equation, applied to the functions ^y 

and ^i , using the values of these functions at si and at s^ . 

The University of Wisconsin 
September, 1916 

dent change of notation would yield another set of formulas precisely similar to those in the 
text. Professor Bliss has pointed out that if an arbitrarily selected set of n- — 1 linearly inde- 
pendent normal solutions is used, formulas can be deduced which are more complicated than 
those in the text, but which apply independently of the $ which is different from zero. 

* Von Escherich, Wiener Berichte, vol. 110 (1901), p. 1376; also Hadamard, loc. cit., 
p. 319, and Bolza, loc. cit., p. 626. 



